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1 Introduction 



The goal of this paper is to show how to produce a piece of rigorous bifurca- 
tion diagram of periodic orbits for an ODE. We study the Rossler system [R] . 
one of the textbook examples of ODEs generating nontrivial dynamics, for the 
parameter range containing two period doubling bifurcations. 

According to the discussion in Kuzniecov textbook |Ku[ Section 2.7] there are 
two extremes in studying bifurcations in dynamical systems. The first one, going 
back to Poincare, is to analyze the appearance (branching) of new invariant 
objects (equilibria or periodic orbits) from the known ones as parameters of the 
system vary. A good reference for this approach is a textbook by Chow and 
Hale |CH) . On the other extreme, it is the approach going back to Andronov 
|An| and Thorn [T] , is to study rearrangements (bifurcations) of the whole phase 
portrait under variations of parameters. It is apparent that the first approach 
is necessarily one of the initial steps in attempting to describe the bifurcations 
in Andronov-Thom sense. In fact in many dimensional systems (even for planar 
maps like the Henon map) achieving the complete description of the phase 
space portrait and its changes appears to be hopeless in view of the results on 
the Henon-like maps [MVl [BCl IWYTl IWY2] . 

While there exists a vast literature on the bifurcation theory, see for exam- 
ple [AAISj ICH| mi IKuj and references given there, and also a lot of numerical 
bifurcation diagrams for various systems can be found in literature (see for 
example references in |Ku] ) . there are virtually no rigorous results on bifurca- 
tions of periodic orbits for ODEs in dimension three or higher in the situation, 
when the periodic orbit undergoing the bifurcation is not given to us analyti- 
cally due to some special symmetries of the system. The basic reason for this 
is: while numerical experiments and/or normal form computations may clearly 
show what is happening (in terms of the bifurcations) we usually lack any rea- 
sonable rigorous estimates about the observed orbits, which prevents us to turn 
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these observations into rigorous statements. To obtain the necessary estimates 
one needs to integrate the variational equations describing the partial deriva- 
tives with respect to initial conditions up to order 3 or higher. This is usually a 
serious problem for rigorous ODE solvers. It turns out that the naive approach: 
applying an ODE solver to the system of variational equations does not work, 
because the methods dealing with the wrapping effect used in the Lohner-type 
algorithms (the most effective rigorous ODE solvers) [Lol IZ11 INJ| break down 
for such system. As the solution of this problem C r -Lohner algorithm has been 
proposed in |WZ| and it is used in the present work. 

Concerning the content of the paper regarding the bifurcation theory itself, 
we were forced to reformulate some well known theorems to make them amenable 
to computer assisted proofs. It is a common feature of all bifurcations theorems 
that the bifurcation point (or rather a candidate) and all necessary data like the 
spectrum and maybe some higher order terms are always given as part of the 
assumptions. But in a nonlinear system we usually do not have explicitly these 
data, in fact the existence of the bifurcation point has to be proved by looking on 
the behavior of the system in some neighborhood. This forces us to reformulate 
some bifurcation theorems in a semi-local way, we have to investigate properties 
of solutions of implicit equations, which are degenerate (due to the presence of 
bifurcations). This is the reason, why from various approaches to bifurcations 
we chose the one developed in |CH] and which is based on the Liapunov-Schmidt 
reduction. 

In our work we focus on the period doubling bifurcation of periodic orbits 
for Rossler equations, in fact we study the Poincare map for Rossler system. 
The paper is organized as follows: in Sections [3l |4] and [5] we discuss the main 
tools used to produce a validated piece of the bifurcation diagram containing 
the period doubling bifurcations. In the remaining sections we give some details 
concerning our results for Rossler system. 

2 Basic definitions 

By N, Z, Q, R, C we denote the set of natural, integer, rational, real and 
complex numbers, respectively. Z_ and Z + are negative and positive integers, 
respectively. By S 1 we will denote a unit circle on the complex plane. 

For K™ we will denote the norm of x by ||a;|| and if the formula for the norm 
is not specified in some context, then it means that one can use any norm there. 
Let x G R s , then B s (x ,r) = {z e M s | ||atr - 4 < r} and B s = B S (0,1). 

Let A : W l — ► M" be a linear map. By Sp(v4) we denote the spectrum of 
A, which is the set of A £ C, such that there exists x ^ C n \ {0}, such that 
Ax = Xx. 

For a map / : X — » Y by dom (/) we will denote the domain of /. For a map 
F : X -v X we will denote the fixed point set by Fix (F, U) = {x E U \ F(x) = 
•'•}• 

Let x = (x\, . . . , x n ) 6 M" . By iti we will denote the projection on i-th 
coordinate, i.e. ni(x) = xi. Analogously for any multiindex a = (ii, 12, ■ ■ ■ , ik) £ 
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Th we define 7r a (x) = (xi lt Xi 2 , ... ,%i k )- Sometimes the points in the phase 
space will have coordinates denoted by different letters, for example z = [y, x, y), 
then we will index the projection by the names of variables, i.e. ivi v ,x){ z ) — x ) 
etc. 

Definition 1 Let f : R" D dom (/) -> R n be C 1 . Let z G dom(/). We say 
that Zo is a hyperbolic fixed point for / iff f(z ) = z and Sp(Df(z ))r\S 1 = 0, 
where D/(zq) is the derivative of f at zq. 

Definition 2 Consider a map f : X D dom (/) — > X. Let x G X. Any 

sequence {x k } keI , where I CZ is a set containing and for any l\ < l 2 < I3 in 
2 if h,h € /, then l 2 G /, such that 

x =x, f(xi) = x i+1 , fori,i + l£l 

will be called an orbit through x. If I = Z_ U {0}, then we will say that {xk}kei 
is a full backward orbit through x. 

Definition 3 Let X be a topological space and let the map f : X D dom (/) — > 
X be continuous. 

Let Z C MP, x Q G Z, Z C dom (/). We define 

W s z {z , f) = {z\ y n >of n (z) G Z, lim P(z) = z } 

n—*oo 

Wz(zoif) = {z I 3 {x n } C Z a full backward orbit through z, such that 
lim x n = z } 

n—> — oo 

W s (z ,f) = {z\ lim f n (z) = z } 

n^oo 

W u (z ,f) — {z\3 {x n } a full backward orbit through z, such that 

lim x n = z } 

n—^—00 

Inv + (ZJ) = {z\V n > f n (z)eZ} 

Inv~ (Z, f) — {z\3 {x n } G Z a full backward orbit through z } 
Inv(ZJ) = Inv + (Z,f)r\Inv-(Z,f) 

If f is known from the context, then we will usually drop it and use W s (zq), 
Wz(zo) etc instead. 

Definition 4 Let P v : M™ — > W 1 , where v belongs to some interval. We say that 
P v has a period doubling bifurcation at (i>o, x ) iff there exists V — [v\, v 2 \ *X C 
IxR", such that the following conditions are satisfied 

• (y ,x ) G mtV, P Vo (x ) = x Q 

• there exists a continuous function Xf p : \v\,v 2 \ — > intX, such that 

Fix{P v ,X) = {x fp ( V )} 
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• there exist two continuous curves Ci : [^0,^2] —> wkX, i = 1,2 , such that 
for v £ [uq, v 2 \ holds 

ciM = c 2 (v ) = Xfp("o) 
ci(v) ^ c 2 {v), v^v a 
P v {cl{v)) = 02(1/), P tf (caM) = ci(i/) 
Fix(P?,X) = { Cl (v),c 2 {v),x fp (v)} 

• i/ie dynamics: 
for v < vq 

Inv{X,P y ) = {x fp {v)} 
For v > vq the maximal invariant set in X Inv(X, P v ) is equal to 

W$(x fp (u),P) n (W^( C1 ^),P^ u w% (c2(v),pS)) 

and is a one- dimensional connected manifold with boundary points c\{v), 
c 2 (»- 

3 Derivation of the conditions for the occur- 
rence of the period doubling bifurcation 

The goal of this section is to present the set of conditions, which guarantee 
the existence of period doubling bifurcation for a given map, and which can be 
verified using rigorous numerics. The main tools used are the Liapunov-Schmidt 
reduction CH and the implicit function theorem. 

Assume that we have a parameter dependent map z 1— > P(y,z), which ap- 
parently undergoes the period doubling bifurcation as the parameter v changes. 
Let Zf p (v) be a fixed point curve for P. We assume that it is regular and we 
can compute it and its all derivatives. 

To prove the existence of the period doubling bifurcation we proceed as in 
|CH| . First we perform the Liapunov-Schmidt reduction to obtain a function 
G:MxMd dom (G) — ► M, whose zeros correspond to fixed points and period 
two points of P v and then we try to describe the solution set for equation 
G(v, x) — 0. Next, through some additional computation of eigenvalues we will 
be able to decide about the hyperbolicity of bifurcating periodic orbits. 

The basic steps of the Liapunov-Schmidt reduction for P 2 are: 

• we choose good coordinates (x, y) £ Rxl"" 1 . It is desirable to choose x in 
the approximate bifurcation direction (in the eigendirection corresponding 
to —1 eigenvalue at the bifurcation point). 



let Z — [ui, v%] x [xi,X2] and Y C M. n 1 be such that the apparent bifur 
cation point (j/q, Xq, yo) belongs to the interior of Z x Y 



4 



• we need to show that there exists a function y(v,x), defined on Z with 
the values in Y, such that 

y-w y (P?(x,y))=0 for (y, x, y) E Z x Y iff y = y(v,x). (1) 

• the bifurcation function G : Z — > R is defined by 

G(i/,x) =1-^^(1, a;))). (2) 

Now, we have to find the solution set of the following equation 

G(v,x)=0, {v,x)eZ. (3) 

Let Xf p (y) = n x (zf p (i>)) be the x-coordinate of the fixed point curve. We 
assume that [1/1,1/2] C dom(x/ p ) and Xf p ([v\, 1/2}) C [xi,X2]. Therefore we have 

G(u,x fp (u)) = 0. (4) 

The idea of solving ^ goes as follows: we introduce a new bifurcation 
function 

9{v,x) = (5) 

X-Xf p {v) 

and then we solve equation g(y,x) = by the implicit function theorem. 

Observe that expression ([5]) defining g{y, x) contains zero in the denominator, 
moreover usually the exact value of Xf p (y) is not known, therefore the formula 
([5]) appears to be useless in rigorous computations. The next lemma will give 
us an integral representation of g, which will not contain any singularities and 
therefore it is well suited for rigorous numerics. 

Lemma 1 Assume F : R™ -» R s is C 1 . Let x, y £ R". Then 

f 1 dF 

F{x)~F(y) = —(t{x-y)+y)dt-{x-y) 

Hence we can define equivalently g : [v\, — > [x\, x%] by 
f 1 dG 

g(v,x) = J —(u,t(x-Xfp(v))+Xf p (u))dt. (6) 

We obtain 

G(v, x) = (x- x f p{v))g(u, x). 
Therefore, we have to determine the solution set of the following equation 

g(v,x) = Q (v,x)eZ, (7) 

where g is defined in 

In the case of the period doubling bifurcation we expect solutions of ([7]) to 
form a regular curve. The following lemma gives a set of conditions, which 
implies this fact. 
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Lemma 2 Let Z = [vi, 1/2] X [a^i,^]- Assume that g : Z — > R is a C k -function, 
k>2. 

Assume that 



d 2 g 








(8) 




> 


0. 




ov 


< 


0, 




(9) 




> 


0. 


for x G [2:1,0:2] 


(10) 


g{v2,xi) 


> 


0. 




(11) 




> 


0. 




(12) 


g(v2,x ) 


< 


0, 


for some xq 6 (xi,X 2 ) 


(13) 



Then there exist Xi,X2, such that X\ < X\ < xq < X2 < X2 and there exists a 
function v : [xi,X2] [v-y, ^2] of class C , such that 

{{v, x) e Z I g(u, x) = } = {(v(x),x),x e [xi, X2]}. 

Moreover, there exists Xq S (xx,x 2 ) such that 

v'{x) > 0, x e (xo, X2) 

v'{x) < 0, x € (xi, Xq) 

v(x) > V\, x £ [xi, X2] 

v{xi) = v(x 2 ) = v 2 . 

Proof: Observe first that from condition ([8|) it follows that for any given v £ 
\y\ , V2] and any c 6 R the equation 

g{v,x) = c, 

has at most two solutions in [xi,X2]. 

From this observation and equations (llllH3| it follows that there exist x~\ 
and X2, such that 

X\ < x~i < Xo < X2 < X2 

{X € [xi,X 2 ] I g{v2,%) = } = {xi,X 2 } 

g{v2,x) > 0, for x < x\ or x > x~2 

g(v 2 ,x) < 0, for x e (xi,x 2 ). 

From the above conditions and conditions © and (fTX))) it follows immediately, 
that there exists function v : \x\,X2\ ^2], such that 

{{v,x) £ Z I g{v,x) = } = {(v(x),x),x € [xi,x 2 ]}. 

By the implicit function theorem function v(x) is of class C k . 

It remains to show the existence of a unique minimum of v{x) and its mono- 
tonicity properties. 

G 



Let y 6 (xi,X2) be any critical point of f(x), i.e v{y) = 0. We will show 
that i)(y) > 0. 

By differentiating twice equation g(y(x), x) = we obtain 

^(v(x),x)(i>(x)) 2 + 2-^jL(v(x), x)i>{x) + 

dg d 2 g 

— (v(x),x)i>{x) + -—?(v{x),x) = 

av ox 1 

Therefore for y we have 

= ^{is(y),y)i>(y) + ^(v(y),y) 

We see that all critical points are strong local minima. This implies that the 
set of critical points consists from just one point. I 

The model for Lemma [2] is given by the function gi (y 1 x) — x 2 — v in the 
neighborhood of point (0,0). By changing signs of v and g we obtain the follow- 
ing model functions gi(y, x) — v + x 2 , g^,{i>,x) = v — x 2 and gi{v, x) = —v — x 2 
for which we can state analogous lemmas. 

Now we can formulate a lemma based on the implicit function theorem 
addressing the assumptions implying intersection of curves solving equation 
G(v, x) = 0, where G arises in through the Liapunov-Schmidt reduction in 
the context of the period doubling bifurcation. 

Lemma 3 Let Z = v%\ x [x±, xs]. Assume that G : Z — > Hi is a C k -function, 
fc>3. 

Assume that there exists a C k -function Xf p : [1/1,1/2] — » (2:1,0:2), such that 
G(v,Xf p (v)) = for v 6 [1/1,1/2]. 
Assume that 

cfiC 

^(Z)>0 (14) 

^(2) + ^(%H^)'[°. 1 ]< ' ( 15 ) 
We assume that following conditions are satisfied for some x\ < Si < 

Xfp(vi) < 62 < X 2 

G(v 1 ,[5 3 ,x 2 ])>0, G{v u [x u 5i]) < (16) 

^MhM) > o, (17) 

G(V2,X2) >0, G{u2,xi) < (18) 
dG 

— ( V2 , Xfp (y 2 )) < (19) 
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Then there exist x\ < x\ < X2 < xi, such that Xf p {v2) G {x\ 1 X2) and a 
function v : [xijica] ~ ¥ [Vi> ^2] of class C , such that 

{(v, x) e Z I G(v, x) = } = C/p u C per = 

{(^,a;/p(i/)), S [2^1,^2]} U {{y(x),x),x G [2Ci,ac 2 ]} (20) 

and i/ie intersection of curves C f p and C per contains exactly one point. 
Moreover, there exists Xq G (x 1,3:2) such that 

v'(x) > 0, x G {xq,xz) 

v'{x) < 0, x G (xi, So) 

1/(1) > ^1, x G [xi, X2] 

v{xx) = v(x 2 ) = v 2 . 

Proof: For the proof we want to apply to Lemma [2] For this end we define g 
as in ([6]). 

We start by showing that (TJJ) and (JTSJ imply that > and |£ (Z) < 

0, respectively. 
We have 



2 S, , f 1 



d 3 G 

— (v, t(x - x fp (u)) + x fp (v))t 2 dt. 



n2 

Hence from ([M)) we obtain immediately that -a-%{Z) > 0. 



^-{v,x)= [ (^-^-{v,t{x-x fp {v))+x fp {v))+ 



dv ' J \dxdv 
d 2 G 



dx 2 



(v, t(x - x fp {v)) + Xf p (v))(l - t)x' fp (v)j dt C 



This and (HSJ) imply that §£(Z) < 0. 

To obtain condition ([TO")) we need to split the interval [xi, X2] into three parts 
[xi, Si], [Si, 82] and [82, x 2 ], so that in the middle part we have the zero of G{yi, ■) 
and we need to use there the integral representation of g. On the remaining 
parts it is enough to verify the signs of G. Hence we see that conditions (fT6T - [l"7|) 
imply ([TO]), 

The remaining assumptions in Lemma follow easily from (p~51 TT^]) . Now we 
use Lemma [21 to obtain function v(x) and condition (|20p . 

It remains to show that curves Cf p and C per defined by (|20[) intersect exactly 
in one point. Observe that these curves intersect because curve Cf p cuts Z into 
two pieces and the end points of the second curve belong to different components, 
which follows directly from the fact that Xf v {i>2) G (x±, X2). 
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Now we turn to the question of the uniqueness of the intersection point. 

Let a, (3 £ [v\, 1*2], a < f3. For t £ [0, 1] let v t = ta + (1 — t)[3 and x t — 
txf p (a) + (1 — t)xfp(P). Observe that for each t £ [0, 1] point (yt,Xt) belongs 
to Z. Let 9 £ {a, [3) be such that x' fp {9) = x fA<*y-*fAP) _ We have 

^(a,x fp {a)) - ^(p,x fp (P)) = 

( d2 ° , w as d2 ° 
— [y t ,xt)(a -p) + - 

\dx^^ tlXt ^ + dtf^^'fp^J dt ) ( a ~i 3 ) c 

Therefore, from above computations and assumption (|15[) it follows that the 
function ^ 1— > T^(v,%fp(v)) is injective on [fi,^]. Observe that from ([S]) it fol- 
lows that, if (f, Xfp{i>)) £ Cf p C~\C per then ^(v, Xf p {v)) = 0, so the intersection 
of Cfp and C per contains at most one point. 1 

Observe that in the above lemma we cannot make the claim that the in- 
tersection point of the curves, which solve equation G(y, x) — is exactly in 
(v(xo),x~o). This can be easily seen in the following example. Let G{v,x) = 
(x — l)(x 2 — v), xi = —2, X2 = 2, v\ = —1 and V2 = 1. It is easy to see 
that all assumptions of Lemma [3] are satisfied, but the intersection of the curves 
{v{x) — x 2 ,x) and (v,x{v) = 1) is not (0,0). On the other hand in the context 
of the period doubling bifurcation the intersection point is (v(xo), Xq), but we 
cannot infer such conclusion from Lemma [3] and we need to use the information 
about the dynamical origin of function G. Now we state the theorem which 
addresses this issue. 

Theorem 4 Let P v : R" D dom(P, y ) — * R™, where v £ I C R be one-parameter 
family of maps of class C k (k > 3), both with respect to the parameter v and 
x £ R n . 

Let Z — [u\, v%] x [x%, X2] and Y C R™ -1 be a closure of open set, such that 
[xx,X2\ x Y C dom(Pj) for v £ \y\,vi). Assume that 

Al for any (i>, x) £ Z there exists a unique y — y(y, x) £ int Y , such that 
y — n y (P 2 (x,y)) = 0. Moreover, we assume that y : Z —>Y is C k . 

A2 there exists C k -function Xf p : [^1,^2] ~^ (a?i , 2^2), such that for v £ [v 1,1/2) 
holds 

Fix(P„,[xi,x 2 ) xY) = {{x fp (is),y( V ,Xfp{iy)))} (21) 

A3 Let 

G{v,x) =x- Tr x (P 2 (x,y(v,x))), for (i/,x) £ Z. 

Assume that G and Xf p satisfy assumptions of Lemma\3[ and let xi,X2, Xo 
and v : [xi,X2] — > [^1,^2] be as in the assertion of Lemma\^ 
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Then the fixed point set of P^ for v G [v\, v 2 \, i.e. 

{(v,x,y)€Z xY | P*{x,y) = (x,v)} 

is equal to the sum of the fixed point set for P v 

Pen = {{v,Xfp{v),y(v,Xf p (v))) \ v e \vi,v 2 \} 

and the period-2 points set 

Per 2 = {(u(x),x,y(v,x)) \ x G [xi,x 2 ]}. 

Sets Peri and Per 2 have exactly one common point {v b , 2%) given by 

(v b , z b ) = (v(x ), (x ,y(i/(xo),x )). 

Moreover, the projections of Peri and Per 2 onto (y, x)-plane have exactly 
one common point (i/f,,a:&) given by 

Proof: From the construction of the bifurcation function G and our assump- 
tions we immediately obtain that 

{(u,x,y) eZ xY I P*(x,y) = (x,y)} = P eri UPer 2 . 

From Lemma [3] we know that projections onto (v, a;)-plane of sets Peri 
and Per 2 intersect exactly in one point, say (v, x). Observe that the point 
[v, x, y{v, x)) belongs to the intersection of Peri and Per 2 . 

It remains to show that (p, x) — (v(xq),Xq). We will show that the function 
x 1 ► v{x) has a local extremum at x. This will imply that x — xq, because by 
Lemma [3] xo is the only local extremum of v(x). 

We reason by contradiction. Let us assume that v'(x) ^ 0. Let U = U v x 
U x x Uy, where U v C [^1,^2], U x C [aj^xa] ana - U y C Y, be neighborhood of 
(v,x,y(v,x)), such that 

P„(x,y) G int ([£1,0:2] x Y), for [y,x,y) G U 

v(a) ^ v(b), for a, b G U x and a ^ b. (22) 

Such U exists because (x,y(P, x)) is a fixed point for Pp and (x,y{p,x)) G 
int ([2:1, £2] x Y). 

Let us take v G U x , such that v ^ x. Then (y, y(y(y),v)) is not a fixed 
point for P v i v \- Points [y, y{y[y),v)) and P v i v \(v, y{y(y),v)) are different, both 
belong to Z and are period-2 points for P v r v )- Therefore they both belong to 
Per2 and 

v{K x P v {v){v,y{v(v),v))) = v(v). (23) 
Observe that from the continuity it follows that 

lim it x P u m (v, y(v(v), v)) = n x P v tg\ (x, y{v{x), x)) = x. 

v^ X 

From the above observation it follows that for v sufficiently close to x points v 
and ic x P v r v \(y, y(u(v), v)) are in U x , but in this situation condition (I23p contra- 
dicts ([22]) . This proves that x = xq. 
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3.1 Hyperbolicity of bifurcating solutions 

The Liapunov-Schmidt projection does not give any direct information about 
the dynamical character of the bifurcating objects. The required information 
concerning the hyperbolicity is of course contained in the spectra of DP V and 
DP 2 and its derivatives. Below we present a lemma addresing this issue. 

Lemma 5 Assume that P u : M™ — > M™ for v G \u\, v 2 \ satisfies all assumptions 
of Theorem^ and in the sequel we will use all the notation introduced there. 
Let z fp (v) = (x fp (v),y{v,x fp {v))). 

fixed points: Assume that there exists e > 0, such that for all v G [i/i, v 2 \ holds 

Sp (z fp (v))\ = A v U B v U 

where X(v) G R has the multiplicity one, A v C {a G C, |a| < 1 — e} and 
B v C {P G C, \/3\ > 1 + e}. Moreover, we assume that 

A(n) c (-1,1) 
A(i*) < -1 

Then the fixed points for P v on curve z f p (y) are hyperbolic for v G \v\ , v 2 \ \ 
{v(xq)} and 

dimW u (z fp (i;-),P l ,-) + 1 = AimW u {z }p { V + ),P u+ ) 

for any v\ < v~~ < v{xq) < v + < v 2 . 

period-2 points: Assume that there exists e > 0, such that on the Per 2 curve 
(i.e. for x G [xi,52]J holds 

I dP 2 \ 
Sp I (x, y(u(x),x)) \=A x yjB x yj { 7 (x)} 

where 7(2:) G M has the multiplicity one, A x C {a G C, \a\ < 1 — e} and 
B x C {[3 G C, |/3| > 1 + e}. Moreover, we assume that 

— (x) < 0, x G [xi,x 2 ] 
0<j(xi) < 1. 

T/ien for x G [2:1, X2] \ {xq} the period two points Zd(x) — (x,y(v(x),x)) 
for P v t x ) are hyperbolic and 

l{%o) = 1 
< j(x) < 1, 
dimW s (z d (x),P^ x) ) = dimW s (z fp (v-),P v -) 

for any v\ < v~ < v{x$) and x G [ah, £2] \ {%o} 
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Proof: The statement about the hyperbolicity of fixed points is obvious. 

For the proof of the second part it is enough to observe that in the bifurcation 
point holds 

\{v(x )) = -l, 7O0) = A(Kzo)) 2 = 1, p-(xo)=0. 



4 Continuation 

To apply the tools described in Section [3] in the part regarding the existence of 
the Liapunov-Schmidt reduction we need to prove the existence and uniqueness 
(locally) of solution of the equation of the form f(a, y) = for a given a, where 
y € R™ and a is a parameter. Similarly, when continuing the fixed point curve 
or period-2 point curve we have solve the existence and the local uniqueness of 
the solution of x — P % {a,x) — 0, where a is the parameter. It turns out that 
both of the above mentioned tasks, can be handled by the same tools. 

In this section we will discuss such tools, the first one consists of classical in- 
terval analysis tools: interval Newton method [A, MolIN] and Krawczyk method 
[Al IKr[ IN] , which can be seen as clever interval versions of the standard Newton 
method. These methods work very efficiently in the situation, where the solution 
sought is well isolated from other solutions and it requires C 1 -estimates, only. 
The second approach, which is based on the implicit function theorem deals 
with situation, when we are close to the bifurcation point and therefore there 
are several solutions close to one another, as in the case of the period doubling 
we have the fixed point and period two points in a small neighborhood. 

4.1 Two methods for proving the existence of zeros for a 
map. 

Let A C R". By [A]i we will denote an interval enclosure of the set A, i.e. the 
smallest set of the form [A]j = [ax, bi] x • • • x [a n , b n ], such that A C [A]i, where 
a h bi eM"U {±00}. 

Theorem 6 (Interval Newton Method JMM^\M) Letx 

C R™ be a convex, 

compact set, f : N — * R" be smooth and fix a point ifAf. Let us denote by 

N(f 1 X,x) = x-[Df(X)]J 1 f(x) (24) 
the Interval Newton Operator for a map f on set X with fixed x 6 X. Then 

• if N(f,X,x) C int X then the map f has unique zero in X. Moreover, if 

is such unique zero of f in X then x* S N(f,X,x). 

• if N(f, X,x) (~\ X = $ then the map f has no zeros in X . 
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Theorem 7 (Interval Krawczyk Method JI\ \Kr\ Utf) Let X C R n be a 

convex, compact set, f : N — ► R n &e smooth and fix a point x £ N . Let C £ R nxrt 
be an isomorphism. Let us denote by 

K{f, C,X,x)=x- Cf(x) + (Id-C- [Df(X)]j)(X - x) (25) 

the Interval Krawczyk Operator for a map f on set X with fixed x £ X and 
matrix C . Then 

• if K(f, C, X, x) C intX then the map f has unique zero in X. Moreover, 
if x* is such unique zero of f in X then x* £ K(f,C,X,x). 

• if K{f, C,X,x) C\ X — % then the map f has no zeros in X . 

4.2 Continuation close to the bifurcation point 

Lemma 8 Assume /„:lx W 1 ^ 1 D X x Y — > W l , v £ [1/1,1/2] is C k function 
both with respect to argument and parameter, with k > 3, such that 

1. for v £ [1/1,1/2] there exists unique fixed point (x/ p (i/),y/p(i/)) for f v in 
X x Y 

2. for all (v,x) £ [1/1,1/2] x X there exists unique y(y,x) £ intT solving 
equation y — Tr y (f„ (x, y)) = and the map y: [1/1, v 2 \ X X —>■ Y is of class 
C k . 

3. the map G{v,x) — x — ir x (f^{x,y(v,x))) satisfies 

d 3 G 

-^(y,x)>0, for v £ [u u i/ 2 ], x £ X (26) 
G(i/,Xi) < 0, G(i/,x 2 ) > 0, forv£\v 1 ,v 2 ] (27) 
Vz/ £ [2/1,1/2] 3x_ £ (xfp(v),x 2 ) G(v,x~) < (28) 

Vi/ £ [1/1,1/2] 3x + E (xi, a;/ p (i/)) G(v,x + ) > (29) 

TTien there exist two C k curves c\ , C2 : [1/1 , 1/2] — ► R n smc/i i/ia£ /or 1/ £ [1/1 , 1/2] 
/ioWs 7Tx(ci(i/)) < Xj p {v) < 7r x (c 2 (i/)) and Cj(i/) is a period two point for f v , 
i = 1,2. 

Moreover, ifforsomev £ {fx, i/ 2 } holds f Uo (ci(v )) = c 2 (i/ ) or f„ (c 2 (vo)) = 
Ci{vo) then for all v £ [v\, v 2 \ 

f v {cM) = c 2 (i/), /,(c 2 (i/)) = ci(i/) (30) 

Proof: The second assumption and |2"6"1) imply that for a fixed i/ the map has 
at most three fixed points in X x Y, From the first assumption we know that /„ 
has unique fixed point (xf p {y), yf P (v)) in IxF. Therefore any zero of G{v, •), 
which is different from (y, x fp{v)) corresponds to a period two point of /„. From 
the continuity of G and from (|27H29| it follows that G has one zero in each of the 
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intervals x per (v) £ (xi,Xf p (u)) and Xp er (u) € (x f p (v) , x 2 ) ■ It is easy to see that 
functions x per are continuous for i = 1,2. We set Ci(y) = (x per (v), y(v, x per {v)). 

We will show the smoothness of x l per , which together with assumption that 
y(x,v) is C k implies the smoothness of c;. It is enough to show that 

dG 

— (v,x per (v))^0, 

because then we can apply the implicit function theorem to obtain the required 
differentiability. Let us fix v S \v\ , v 2 \ ■ Observe that condition ([26]) implies that 
for any fixed v the function x 1— » ^-(y, x) has at most two zeros in [x\, X2]. From 
remaining assumptions it is clear that on interval [x\, Xf p (i/)] and [xf p (v),x 2 ] 
the function x 1— > G(v, x) has strictly positive maximum and strictly negative 
minimum, respectively. Therefore these extremal points are zeros of ^(v,x) 
and obviously they are different from points x per (v\ which are zeros of G{v 1 ■). 
Hence we have shown that ^:{v,x per {y)) ^ 0. 

Assume that vq = v\ (the other case is analogous). From the implicit func- 
tion theorem it follows that for some v 1 > v\ condition (f50| is satisfied for 
v\ < v < v' . Let v m be supremum of such v 1 < It is easy to see that 
v m = v%i because at v m is also satisfied by the continuity and implicit function 
theorem allows us to extend the range of v satisfying (|30|) to the right if v m < v 2 - 



5 Extracting the dynamical information from Lia- 
punov- Schmidt reduction 

As was mentioned already in Section [3. II the Liapunov-Schmidt projection does 
not give us any direct information about the dynamics of bifurcating solutions 
regarding the invariant manifolds of the bifurcating objects as required by Dcf.|H 
In this section following the ideas of de Oliveira and Hale [Hj I OH] , we show 
that the information obtained from the Liapunov-Schmidt reduction and the 
spectrum of the bifurcating fixed point curve is enough to say precisely, what is 
the dynamics in the neighbourhood of the bifurcation point. 

Our argument follow the ideas from |CH[ Chapter 9, Thm. 3.1 and 4.2], 
where an analogous problem was considered for fixed points for ODEs and pe- 
riodic orbits for periodically forced ODEs. The notion of the central manifold 
[K] plays crucial role in this proof. 

Theorem 9 Let P u : W n -> E n for v e [vi,v 2 ] be a C k -map (k > 3) both with 
respect to v and its arguments. Assume that on the set V = [v\, v 2 \ x [£1,3:2] xY, 
where Y C R" _1 is a closure of an open set, we were able to perform the 
Liapunov-Schmidt reduction and verify assumptions of Theorem^ Let (fb,Zb) 
be the bifurcation point and Zf p (y) = (xf p (v),y(i>,Xf p (i>)) be the fixed point 
curve for P v in V . 
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Let v be the eigenvector of (zb) corresponding to the eigenvalue — 1. We 
assume that ir x v 0. 

Assume that there exists e > 0, such that for all v G \y\, 1/2] holds 



where \{v) G R has the multiplicity one, A v C {a G C, |a| < 1 — e} awe? 
-B„ C {f3 G C, |/3| > 1 + e}. Moreover, we assume that 



Then the map P has a period doubling bifurcation at (i/f,, x;,, z/(z^,, Xf,)). 

Proof: Let v : [xi,^] ~~ > ^2] be the function from assumption A3 of Theo- 
rem [4] (in fact of Lemma [3]) is satisfied. In the notation used in Theorem [4] we 
have (i>b,Zb) — (z/(xo), (xo, y(y(xa), So)))- Let c\(y) and C2(y) be respectively 
lower and upper branch of the graph of the function x — > f(x) giving period-2 
points - see Fig. [1] 

Let us define a map : V — > R x R" by 



Consider the spectrum of DH(vb, Zb). It is easy to see that +1 is an eigenvalue 
of DH(vb, Zb) of multiplicity one, A = — 1 has also multiplicity one and all other 
eigenvalues are off the unit circle. 

We apply the center manifold theorem [Kj [G] IHPSj to H in the neighbour- 
hood of {vb,Zb). Therefore, there exists a neighbourhood M of (f&,Zfe) and 
two-dimensional center manifold W c C M such that 



W c is tangent at {i>b,Zb) to the subspace spanned by vectors {(1,0), (0, v)} C 
R x 1". Observe that from our assumption about v, i.e. ir x (v) / 0, it follows 
that we can use on W c in the neighbourhood of (vb, Zb) the same coordinates 
{v, x) as in the Liapunov-Schmidt reduction. There exists a neighbourhood 
of (vb,Zb) denoted by U = [^1,^2] x [xi,X2] x Y C M Pi V and C fc -functions 
h : [vi, V2] x [xi, X2] — > Y and / : [z/i, vi\ x [xi, X2] — > R satisfying 



Let us stress that the dynamics of P v in W c is one-dimensional , namely 
that of x 1— > /(^, x). 




AK) c (-1,1) 
A(f 2 ) < -1. 



H{u,z) = {u,P{u,z)). 



V(v, z) G W c if z) G M then #(1/, z) G W C) for i = -1, 1 

Inv (M, #) C I¥ c 



VF C = {(^x,%,x))} 
P(v,x,h{v,x)) = {f(v,x),h(v,f(v,x))) 
Inv ({/, ff) C W c . 



(31) 
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Figure 1: Location of sets Ai, Bi with respect to zeros of G. 



From the Liapunov-Schmidt reduction we know that a point (v,x,y) G U 
has period one or two with respect to map H iff y — y(v, x) and Giy, x) = 0. Let 
N = [i/i, v 2 \ x [xi, x 2 ]. If U is chosen to be sufficiently close to the bifurcation 
point, then the set N\{(v,x) | G{v, x) = 0} has four connected components - 
see FigQ] Namely, 

Ai = {(v,x) e N | {{v < Ub) and (a; < Xf p (v))) or 

((u > v h ) and (a; < ci(v)))} , 
Ai = {{v,x) £ N | {{v < Ub) and (a; > Xf p {u))) or 

((u > v h ) and (a; > c 2 (f)))} , 
Bi = {{v, x) e AT | (f > and (xf p (v) > x > ci(v))} , 
S 2 = {(f, x) e N | (i/ > v b ) and (x/p(v) < x < c 2 (v))} . 

We also require that 



On each of these components the function d{v,x) — x — f(y,f(y,x)) must 
have a constant sign. Observe that on A 2 we have 



because Zf p (Pi) is attracting on W c and we consider the second iterate. Anal- 
ogously we obtain 



xi < ci(i> 2 ) < c 2 (i> 2 ) < i 2 - 



(32) 



x- f{v,f{v,x)) > 0, for (u, x) e A 2 



(33) 



f{v,f(v,x)) < 0, for (z/,x) e Ai 



(34) 



For the component B 2 we have 



x ~ f{v,f{v,x)) < 0, for 0, x) e B 2 
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because Xf p (u 2 ) is repelling on W c and we consider the second iterate. Analo- 
gously 

x-f(vj(v,x)) >0, for (v, x) £ B t . 

For a subset Z C N by Z„ we will denote Z v = \x : (y, x) £ N}. Observe 
that for each v £ [Di, P2] holds 

/*(0A*)„)n[zi,£ 2 ] C (A,),, fl{{Bi) u ) n [xi,i 2 ] C i = 1,2. (35) 

For the proof of (|35p observe that map l(v,x) — (u, f(v,x)) maps connected 
components of N \ {G(v, x) — 0} into connected components, i.e. for any 
5" G {A 1 ,A 2 ,B 1 ,B 2 } there exists T = T(S) £ {A ll A 2l B 1 ,B 2 }, such that 

(l(S)nN) CT(S), (36) 

because 

KG^iQ) f)N) n N = g _1 (o) n N = r 1 (G -1 (0) DN)nN 

Observe that the relevant eigenvalue of DP V at zj p (v) describing the dynamics 
on W c is A(f ), which is real and since we consider the second iterate we see that 
in the neighbourhood of the fixed point curve we have points mapped into the 
same component. This together with (I36[) proves (|35p. 
From the above considerations we obtain for v < vi, 

x fp( v ) < /(^J x )) < x > for x fp( v ) < x < x 2 
x fp{ v ) > /(^ f{ v , x )) > x i f° r x fp v ) > x>x\. 

The above conditions, (|3ip and nonexistence of other period two points in U 
imply that 

Inv (ir Xj yU, P y ) = {z fp {v)}, for v £ [v, v b ]. 
For v £ (yi,, v 2 \ we have 

x fp( v ) < x < f( v i /(^ x )) < c 2{v), for Xf p (v) < x < c 2 (v) 
x fpi v ) > x > f{v,f{v,x)) > c\{v), for Xf p {v) > x > c\{v) 

< f(v, f(v,x)) < x, iorx>c 2 (v) 
ci(f) > f(v, f(v,x)) > x, iovx<ci{v). 

The above conditions, conditions (|32l3ip and nonexistence of other period two 
points in U imply that for v £ [y^, v 2 \ 

Inv (Tr X} yU, P u ) = {(x,h(v,x)) \ x £ [a (v) , c 2 (v)] } 



We would like to stress here that, contrary to all previous theorems and 
lemmas, in the proof of the above theorem we prove the statements about the 
invariant manifold of bifurcating orbits from Definition [4] on some set U, whose 



17 



size we do not control, whereas it is given by the range of the existence of 
the central manifold. In principle, this range can be inferred from the proof 
of the center manifold theorem, but it will be an interesting task to develop 
a computable approach, which will allow to rigorously prove these facts on 
the whole set V. Such task will require explicite estimates about the central 
manifold in the region very close to the bifurcation and some other tools, may 
be of Conley index type [MM| , further away from the bifurcation. 

6 Application to the Rossler system. 

Consider an autonomous ODE in K 3 called the Rossler system [R] 



The classical parameter values (considered by Rossler) are a = 5.7 and b = 0.2. 
For the remainder of this paper we fix b = 0.2. 

The system (j37]) has been extensively studied in the literature numerically 
and is treated in the literature as one of classical examples of systems generating 
chaotic attractor. Yet, the number of rigorous results concerning it is very small. 
In Fig. [2] we show a numerically obtained bifurcation diagram for periodic orbits 
on section x — with b = 0.2 and a as parameter. We see that when the 
parameter a increases from 2 to 5.7 one observes a cascade of period doubling 
bifurcations. In Fig. [3] we show some periodic orbits for different values of a. 
Our goal in this section is to validate the part of the bifurcation diagram in 
Fig[2] containing two first period doublings using the approach introduced in 
the previous sections. 

Let us list the few known rigorous results about (|3"T|) . Pilarczyk (see [P] and 
references given there) gave a computer assisted proof of the following facts: for 
a = 2.2 there exists periodic orbit, for a = 3.1 there exists two periodic orbits. 
However from his proof one cannot infer any information about the dynamical 
character of these orbits. He constructs suitable isolating neighborhoods, which 
have an index of an attracting or a hyperbolic orbit with one unstable direction, 
but no such claim can be made about the periodic orbit proved to exists. In 
fact we do not even known, wether this orbit is unique. 

Finally, for the classical parameter values b — 0.2 and a — 5.7 the system is 
chaotic |Z3j in the following sense: a suitable Poincare map has an invariant set 
S and the dynamics on S contains the shift map dynamics on three symbols. 

Before proceeding any further we need to introduce some notation. Let 
n = {(x,y,z) £ M 3 | x = 0, x 1 > 0} be a Poincare section. Since for u E II 
the first coordinate is equal to zero we will use the remaining two coordinates 
u = (y, z) to represent a point on n. For a fixed parameter value a > by P a = 
{P a , y , Pa,z) '■ n — > n we will denote the corresponding Poincare return map. By 
P we will denote the map defined by P(a,y,z) = (a, P a ,y(y, z), P a)Z (y, z)). 




(37) 
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3.0 3.5 4.0 4.5 5.0 5.5 

Figure 2: Bifurcation diagram for the Rossler system 




-4 -2 2 4 -4 -2 2 4 




Figure 3: Periodic orbits corresponding to fixed point, period two point, period 
four point and period eight point for the Poincare map. Parameter values are 
a = 2.8, a = 3.5, a — 4 and a = 4.2. 



19 



Apparently the first period doubling bifurcation is observed for a « 2.832445 
and the second one for a « 3.837358. In the remainder of this section we discuss 
the computer assisted proof of the existence of both these bifurcations. In our 
presentation we will discuss the first one more in details, while for the second 
one we will just state relevant lemmas and estimates. 

Let uo — (yo, zq) be an approximate fixed point for P a „, i.e. we set 

mo = (2/0, zo) = (-4.7946653021070986256, 0.052488098609082899093) (38) 
and put 

I" 0.99999765967819775891 -0.9582095926217468751" 
- [0.0021634782474835700244 -0.28606708410382636343 

The columns of M are normalized approximate eigenvectors of DP aa (uq), where 
first column corresponds to the eigenvalue close to —1 and the second one to 
the eigenvalue close to zero. On section II we choose new coordinates (y, z) = 
M _1 ((t/, z) — u ) and since, in the sequel, we will use only the new coordinates 
we will drop the tilde. 
Define 

A = [oi,a 2 ] = [2.83244,2.832446] 

Y = [yiM = 1-3107 -[-1,1]-10- 3 
Z = [z!,z 2 ] = 1.3107- [-1,1] • 10~ 4 

Now our goal is present the proof of the following theorem 

Theorem 10 The map P a has a period doubling bifurcation at some point 
(a,y, z) e int (A x Y x Z). 

Remark 11 The existence of period doubling bifurcation is a local phenomenon. 
In fact the sets A, Y , Z can be chosen to be smaller which speed up the proof 
(13 minutes versus 87 minutes), namely we were able to prove the existence of 
period doubling bifurcation in the set 

A = [a u a 2 ] = [2.83244,2.832445028] 

Y = [ yi ,y 2 ] = [-l,l]-10- 4 
Z = [ Zl ,z 2 ] = [-1,1] -10- 5 

However, the choice of larger set facilitates the proof of the existence of con- 
necting branch of period two points between first and second period doubling 
bifurcation, because decreasing a 2 results in the eigenvalue of period-two points 
to be very close to 1, which makes it very difficult to rigorously continue it. 

6.1 The existence of fixed point curve. 

Lemma 12 There exists function (yf p , Zf p ) : A — > Y x Z of class C°° such that 
for (a,y, z) G A x Y x Z holds 

Pa{y, z) = (y, z) iff (y, z) = (y fp (a),z fp (a)) 



(39) 



20 



and 

y' fp (A) c [-1.3336825610133946629, -1.3275439332565022177] (40) 

Proof: The proof, which is computer assisted, consists from two parts, in the 
first one we prove the existence of the fixed point curve and in the second part 
we establish estimate (14U1) . 

For the first part, we use the Interval Newton Method (Theorem [5]) and 
C^-Lohner algorithm to prove that for a G A there exists a unique fixed point 
(yfp(a), Zf p (a)) for P a in Y x Z. In computations we insert the whole set AxYx 
Z as an initial condition in our routine, which computes the Interval Newton 
Operator and obtain that the for all a G A the fixed point (yf p (a),Zf p (a)) 
belongs to the set 



N :=N(ld -P a ,Y x Z) = 

-2.838378938597049559,3.2727784971172813446] ■ 10~ 5 ' 
-4.8121450471307824034,4.2979575521536656939] • 10~ 6 

To obtain (|4D|) we apply C 1 -Lohner algorithm |Zlj to the system 

x' = —y — z 
y' = x + by 
z' = b + z[x — a) 
a' = 

with b = 0.2 in order to compute a bound for y'j (A). Differentiating 

P(a,yjp(a),z fp (a)) = (a,y fp (a), z fp {a)) 
with respect to a we obtain 

d Pa,,y ( 1 _ dP a , z \ dPa.z 0P a , y 



\ dy 



dz J dz dy 



(41) 



(42) 



/ iyu. y dz J da dz 

Vfp = fnp~ ■ A f dP a , z ^ dP a , y dP a - ' 4 ' ! ! 



where the partial derivatives of P are evaluated at (y/ p (a), z/ p (a)). 

We use the set Ax N, where N is defined in (14111 as initial condition in our 
routine which computes partial derivatives of P and after substituting them to 
(|4*3")) we obtain a bound for y'j p as in (|4*0"]). 

We used the Taylor method of order 14 and the time step equal to 0.02 to 
integrate the system (|37| in R 3 for the first part of the proof and the order 
10 and the time step 0.01 when we integrate the extended system (j42|) in the 
second part. I 
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Lemma 13 The eigenvalues Ai,A2 : A — > K of DP a (yf p (a), Zf p (a)) are given 
by 

1 fdP a , y , dP a , 



1 ^P a , y , 9P a , 



A 2 (a) = -(-^ + -^ + .(«)], 



7&p„,„ 5P Q ,A 2 A dP a , y dP ar 



•to = iR 4 -^ +4 



\ dy dz J dz dy 

where partial derivatives of P are evaluated at (yj p (a),Zf p (a)). Let v(a) be the 
normalized eigenvector corresponding to eigenvalue Ai(a). Then 

Ai(ai) e [-0.99999781944914578613,-0.99999548919217751131] 

Ai(a 2 ) 6 [-1.00000064581599335217,-1.00000064581598072628] 

A 2 (A) c [-0.0013533261367103342071,0.0013530378340487671934] 

A; (A) c [-0.70107900728585614836,-0.62770519734197127715] 

i>„(A) C ±[0.99728887963031764841,1.0027184248801992439] 

where v y denotes the y coordinate of v. 

Proof: We leave the derivation of formulas for Ai, A2 to the reader. We used the 
C^-Lohner algorithm applied to the system (|37|) in order to compute bounds for 
Ai(ai) and Ai(a2). Since the parameter 02 has been chosen to be very close to 
the bifurcation parameter we find difficulties with the verification of condition 
Ai(ci2) < — 1 in computations performed in interval arithmetics based on dou- 
ble precision (52-bit mantissa) boundary value type. In our computations we 
used interval arithmetics based on float numbers with 150-bit mantissa (MPFR 
|MPFR| and GMP [GMP] packages). 

Since the eigenvalue Ai(a) of DP a (yf p (a), z/ p (a)) is given by an explicit for- 
mula one can express A^ (a) in terms of first and second order partial derivatives 
of P. We obtain 

, s \ ( d dP a , y d dP a:Z , 

Al(a) = 2UlT + " " s(fl) 



s'(a) = 



1 fdP a ,y dP a A ( d dP a ,y d 8P a , 



s(a) \ dy dz J \da dy da dz 

2 ( ( d dP a ,y \ dP a , Z ( d dP a A dP a ,y 



s(a) \ \da dz ) dy \da dy J dz 
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order 


step 


X 2 (A), X[(A) 


10 


0.03 


Ai(ai) 


10 


0.1 


^1(02) - 150-bit precision 


14 


0.05 



Table 1: Parameters of the C 1 — C 2 -Lohner algorithms. 



where the symbols -§^—q^ 
d 0P a 



— dPa ' y and -§- 9P ,"'* should be understood as 



d 2 P a . 



da dy 



iy f p(a),z fp (a)) 



dady 

9 2 P a z 

dy 2 (yfp( a )' z fp( a ))y'f P ( a ) 

d 2 P a z 

-g-^-(Vf P {a),z fp (a))z' fp (a) 



da dy 

(y fp (a) : z fp (a)) = 
+ 
+ 



and y'f p (a) and z' fj> (a) can be computed as in (|4"3"1) . Next, we applied the C 2 - 
Lohner algorithm |WZ] to the extended system (|42f in order to compute a bound 
for the first and the second order partial derivatives of P and in consequence a 
bound for X'^A). 

We inserted A x N, where N is defined in (|41|) . as the initial condition 
in our routine, which computes the partial derivatives of Poincare map up to 
second order. In these computations we simultaneously computed bounds for 
X2(A) and X[(A). The parameter settings of the Taylor method used in the 
computations are listed in Table [1] 



6.2 The existence of Liapunov-Schmidt reduction. 
Lemma 14 For all {a,y) £ A x Y there exists unique z = z(a, y) £ Z such that 

Pa, z (V, z)=z iff z = y) (44) 

and the map z: A x Y — > Z is smooth of class C°° . Moreover, the map G : A x 
Y — > R defined by 



satisfies 

d 3 G 
dy 3 



G(a, y)=y- Pl y (y, z{a, y)) 



(AxY)c [1.8296823158090675943, 7.2204769494502958338] (45) 



d 2 G 

dy 2 



(A x Y) c [-0.2084557586786322414, 0.2080871792788867581] (46) 
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Proof: Let us fix (a,y) E A x Y and define a function V a , y : Z — » K by 
T/ Qjy (z) = z ~ P 2 z (y, z). The computer assisted proof of this Lemma consists of 
the following steps 

• We divide interval A onto 30 parts. For each subinterval A in this covering 
we proceed as follows 

• Using Interval Newton Method (Theorem[6]) we verified that for all (a, y) G 
A x Y the function V a , y has exactly one zero in Z . Denote this zero by 
z(a,y). This defines the unique map z: A x Y ^ Z which is smooth by 
implicit function theorem and which satisfies (|4"4")l . 

• Let Z denote a bound for z(A, Y) resulting from the previous step. Dif- 
ferentiating z(a, y) — P 2 z {y, z{a,y)) = with respect to y we obtain 

dPg, z 

dy 

d 2 Pg, Z dZ a + jPP a , z dz + d 2 Pg, Z 

dz 2 dy dydz dy dy 2 

d 3 Pg, z (dzV | z d 3 P a , z fdzV 

dz 3 \dy J dydz 2 \dy J 

z fd 2 z d 2 p a , z | d 3 p a . z \ Ch 

\dy 2 dz 2 dy 2 dz J dy 

z d 2 z d 2 Pg, Z | d 3 P a>Z 

dy 2 dydz dy 3 

We see that we can compute all the partial derivatives of z(a, y) as a 
functions of partial derivatives of P. Hence, partial derivatives of G(a, y) — 
V ~ Pa,y(Vi z i a i y)) can be expressed in terms of partial derivatives of P. 
Using the C 3 -Lohner algorithm |WZj applied to the system (|37|) with a 
range of parameter values A and an initial condition Y x Z we computed 
bounds of partial derivatives of Poincare map P up to third order and an 
estimation for f^r(A x Y) and ^p-(A x Y). The estimates (gSJ and p6t 
are an interval enclosures of the estimates obtained in each of 30 steps. 

We used 6-th order Taylor method with the time step 0.04, both, to verify 
the existence of z(a, y) and to compute higher order partial derivatives of P. I 



dPg, z \ ^Z 

dz J dy 
dP a A d^ 
dz J dy 2 

dPg, z \ ^ 

dz J dy 3 
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order 


step 


grid 


remarks 


^(a 2 ,y fp (a 2 )) 


14 


0.05 




150-bit mantissa 


G{a 2 ,yi) 


14 


0.05 




150-bit mantissa 


G{a 2 ,y 2 ) 


14 


0.05 




150-bit mantissa 




6 


0.05 


5 x 30 


integration of (H2")) 


f(W}xy) 


10 


0.05 


1 x 16000 


noncqual parts 



Table 2: Parameters of the C° — C 2 -Lohner algorithms. 



6.3 The existence of period doubling bifurcation for P. 
Lemma 15 For A = [01,02], Y = [y\, y 2 ] the following estimations hold true 

[-1.2916325, -1.2916323] • 10~ 6 (47) 

[-1.15, -1.07] -lO" 13 (48) 
[5.2,5.21] ■ 10" 12 (49) 

[-2.421398492231531, -0.278863623843693] (50) 

[0.83,16.87] • 10" 6 (51) 

Proof: The estimations have been obtained using C° — C 1 — C 2 -Lohner al- 
gorithms applied to the systems ([37]) and (|42|) . The verification of conditions 
(|47H49|) required computations in interval arithmetics based on 150-bit mantissa 
floating points. 

The settings of C°-C 2 -Lohner methods for the above computations are listed 
in Tabled I 

Proof of Theorem HOt The assertion follows from Theorems IH [9] and 
numerical Lemmas fTSl H3 HH IT51 

Indeed, assumptions of Theorem [4] has been verified in 

• Al - Lemma [T4l 

• A2 - from Lemma 1121 there exists a fixed point curve (yf p ,Zf p ): A — > 
(2/1 j 2/2) and from Lemma fT4l it has form as desired in A2 

• A3 - i 7$ (A x Y) because of g5J. 

From j46l), (HDD and gO]) it follows that i ^(AxY) + 0y' fp (A)-[O, 1]. 

Finally, Lemma[l5]guarantees that the remaining assumptions of Lemma[3] 
with £1 = e v = +1 and [81,62] = [j/ijjfe]- 



dy 



(a 2 ,y fp (a 2 )) G 



G{a 2 ,yi) g 
G(a 2 ,y 2 ) G 



d 2 G 
dyda 



(A x y) c 



^({«i}xr) c 
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Finally, from Lemma [T3] we see that the assumptions about the spectrum of 
DP a (A) and an eigenvector v(a) as desired in Theorem [9] are satisfied. I 

6.4 The existence of second period doubling bifurcation. 

In Section [63] we gave a computer assisted proof that for some parameter value 
a\ G [2.83244,2.832446] period doubling bifurcation occurs for P &1 - In this 
section we use similar arguments in order to prove that Pg has period doubling 
bifurcation for some d 2 G [3.83735812, 3.837358168411]. 

Since the arguments used to prove the existence of second period doubling 
bifurcation are the same as in the first period doubling bifurcation we omit the 
details and we present only the sets and the necessary estimates. 

Define 

A 2 = [03,04] = [3.83735812,3.837358168411] 
Y 2 = [2/3,2/4] = [-1.1, 1.1]- 1(T 6 
Z 2 = 3^ 

u 2 = (-4.5003284169596655673, 0.043136987520848421584) 

I" 0.99999908059259889903 0.82277742767392003653 " 
2 ~ |_0.0013560287448822982113 -0.56836370794614177182 

The point u 2 is an approximate period two point for parameter value 04, 
and the columns of matrix M 2 are normalized eigenvectors of DP% , where the 
first column corresponds to eigenvalue close to —1. 

On the Poincare section LT we will use a coordinates (y, z) = M 2 _1 (-u — u 2 ), 
where u denotes a point in cartesian coordinates. In this subsection we will use 
only these coordinates. 

Theorem 16 The Poincare map P 2 has a period doubling bifurcation at some 
point (a 2 , 2/2, 22) G int (A 2 x Y 2 x Z 2 ). 

The proof is a consequence of the following lemmas (proved with computer 
assistance) 

Lemma 17 There exist function (y per , z per ) : A 2 — > Y 2 x Z 2 smooth of class 
C°° such that for (a, y, z) G A 2 x Y 2 x Z 2 holds 

p a(y, z ) = [V, z ) iff (V, z ) = (y P er(a), z per (a)) 

and 

y'(A 2 ) c [-0.36435039423614490328,-0.36419313389173590956] 



26 



Lemma 18 Let Ai, A2 : A — > R be eigenvalues of DP 2 (y per (a), z per (a)) defined 
by similar formulas as in Lemma \1S\ Let v(a) be the normalized eigenvector 
corresponding to eigenvalue Ai(et). Then 

Ai(a 3 ) e [-0.99999992011934590863,-0.99999992005484927837] 

Ai(a 4 ) £ [-1.00000000000149573159618, -1.00000000000149573159615] 

MA 2 ) C [-7.7304566166653588839, 7.7302177026359675856] • 10" 5 

X[(A 2 ) c [-1.6554066232416912996,-1.6460891324715511974] 

TTyv{A 2 ) c ±[0.99984657385734598822,1.0001534381577519284]. 

Lemma 19 For all (o, y) G A 2 X Y 2 there exists unigue z = z(a,y) £ Z 2 such 
that 

p a,z(y> z ) = z iff z = z{a,y) 
and the map z: A 2 x Y 2 — > Z 2 is smooth of class C°° . Moreover, the map 
G : A 2 x Y 2 -> R denned 6y 

G(a, y) = y- Pa, v (y, z(a, y)) 

satisfies 

d 3 G, 



dy 3 
d 2 G 



(A x Y) c [11.780861336872181511,22.544626008881969881] 
(A x Y) C [-0.12474597648618415136,0.12474408945310766494] 



dy 2 

Lemma 20 The following estimations hold true 
dG 

— (a 4 ,y per (a 4 )) € [-2.992, -2.991] • 10" 12 

G(a 4 ,y 3 ) e [-5.13,-5.11] ■ 10" 19 

G{a 4 ,y 4 ) € [5.21,5.22] • 10" 19 

a 2 G 



dyda 



(A 2 x F 2 ) C [-4.354355790892265432, -2.244876361570084633] 



9G 

— ({a 3 } x F 2 ) c [0.99, 30.98] • 10~ 8 
dy 

Parameter settings of computations involved in proofs of the above lemmas 
are listed in Table G3 



7 Continuation of bifurcation diagram 

In the previous sections we proved that the map P a has period doubling bifur- 
cations for parameter values a\ 6 [01,02] = A and a 2 £ [03,04] = A 2 in sets 
Y x Z and Y 2 x Z 2 , respectively. 

Our goal now is to connect these bifurcations with the curve of period two 
points. More precisely, we prove the following result 
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order 


step 


grid 


remarks 




10 


0.04 


150 x 1 






10 


0.04 


150 x 1 




-^-(a4,y P er(a4j) 


25 


0.05 




150-bit mantissa 


G(a 4 ,y 3 ),G(a 4 ,y4) 


25 


0.05 




150-bit mantissa 




10 


0.05 


150 x 1 


integration of (|42|) 


f±({« 3 } x y 9 ) 


14 


0.05 


1 x 10000 





Table 3: Parameters of the C° — C 3 -Lohner algorithms in the proof of the 
existence of second period doubling bifurcation. 

Theorem 21 There exists a continuous curve 

(Vper, Zper) (fll, 62] M 2 

of period two points for P a . Moreover, 



(y per (a), z per (a)), P a {y P er{a) 7 z per (a)) G Y x Z for a x < a < a 2 
(yper(a),Zp er (a)), P%(y per (a), z per (a)) G Y 2 x Z 2 for a 3 < a < a 2 . 

Therefore curve (y pe r, z per ) connects the two bifurcation points for a = a\ and 
a = a 2 . 

The proof of the existence of a branch of period two points for P consists of 
the following steps. 

1. the existence of continuous curve of period two points on intervals (ai, a 2 ] 
and [03, a 2 ] is a consequence of Theorem [TU] and Theorem ll61 respectively. 

2. for parameter values slightly above a 2 , a 2 < a < a, with a — a 2 small, 
we extend this curve using Lemma [HI which requires some C 3 estimates 
(hence it is demanding computationally). 

3. for parameters far from a 2 up to (13, i.e. a < a < 03, we verify the existence 
of period two point curves using Krawczyk method (Theorem [7]), which 
requires only C 1 estimates. 

4. Since we use different methods for proving the existence of segments of 
period two points curve over some intervals in [01,02] it is necessary to 
verify that these segments can be glued to produce continuous curve. 
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At first, it appears that step 2, requiring costly C 3 computations, is not neces- 
sary, because in step 3 we can consider also points close to a 2 using C 1 com- 
putations. But it turns out that, while in principle possible, this approach may 
require very large computation times, because the hyperbolicity is very weak 
there, due to the fact that one eigenvalue of P% 2 is very close to 1. 

To deal with this problem we used Lemma [5] to prove that for parameter 
values slightly above ai there exists a continuous branch of period two points. 
Algorithm [1] is designed to verify assumptions of Lemma [5J In Lemma [22] we 
prove its correctness. 

Definition 5 Let U C R n be a bounded set. We say that Q C 2 R ™ is a grid of 
U if 

1. Q is a finite set and each G G Q is a closed set 

In our algorithms, which will be presented below, we always use grids con- 
sisting of interval sets, i.e. sets which are cartesian products of intervals, most 
of the time uniform grids, which are defined as follows. 

Definition 6 LetY = H™ =1 Yi, where Yi — [ai,bi\ forai < bi and let (g%, . . . ,g n ) G 
Z+. We define a (uniform) giXg 2 X- ■ ■ xg„-grid for Y denoted byQ{gl t . . . ,g n ,Y) 
as follows. 

For any (jx, . . . ,j n ) G Z +; such that ji < g^ we set 



O'i , ai -r Ji 



9ji,...,j n — 

Then G(g^, • ■ ■ , g n , Y) is a collection of all gj 1 ,...,j„ ■ 



(52) 



Lemma 22 If Algorithm [I] is called with its arguments v\, v%, g±, g2, <?3, g x , t, 
X x Y and f v and it does not throw an exception then the assumptions of 
Lemma\$ o-re satisfied for f v ,v^\u\^v^\ on X x Y . 

Proof: The assumption about existence of fixed point curve is verified in lines 
15-19 since for all v G [fx, ^2] the Interval Newton Operator satisfies assump- 
tions op Theorem El 

The existence of Liapunov-Schmidt reduction together with condition (I26p 
is verified in lines 2-8. In lines 4-6 we see that y(v,x) which solves equation 
V ~ Ky{fv( x i v)) i s unique for fixed z/, therefore by the implicit function theorem 
y(y , x) is smooth and we can compute map G and its partial derivatives. 

In lines 9-10 we verify that G(f2, min(X)) < and G(f2, max(X)) > 0. 
Since |f (v, min(X)) > and (f, maxpQ) < (lines 11-14) we see that for 
v G [ui,u 2 ] holds G{v,vam(X)) < and G(f,max(X)) > 0. Therefore IpTj) 
holds true. 
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Algorithm 1: verification of assumptions of Lemma [8] 

Data: [vi, v 2 \ - an interval, gi, g 2 , g%, g x - integers, t 6 (0, 1) - float 

number, X x Y - a convex, compact set, f v - parameterized family 
of maps 

Result: If algorithms stops and does not throw an exception then 
assumptions of Lemma |8] are satisfied 

begin 

Q\ < — .9i x %-grid for v 2 ] X X; 
foreach v x X S Q\ do 

IntervalNewtonOperator(Id y — iry ° f£{X, •), y, center(F)); 
if not C inty then 

I throw Liapunov-Schmidt reduction not verified 

if not ^§-(p,X) > then 
|^ throw condition (|26| is not satisfied 

if (not G(i/ 2 ,min(I)) < 0) or (not G{v 2l max(X)) > 0) then 
I throw condition ([27| is not satisfied 

02 < — .92-grid for [v x , v 2 ]; 
foreach v £ Q 2 do 

if (hoi ^(P,a;i) > 0j or (not ^(v,x 2 ) < 0) then 
|_ throw condition (|27|) is not satisfied 

03 < — .93-grid for v 2 ]\ 
foreach v 6 £3 do 

(X, Y - ) < — IntervalNewtonOperator(Id -^,1x7, center(X x7)); 
if not (X, Y) C int {X x y) then 
I throw fixed points curve not verified 

x + < — (1 - t) min(X) + t min(X); 
X- < — (1 — t) max(X) + t max(X); 

if (not G(min(p), x+) > 0) or (not G(min(£ i ), X-) < 0) then 
|_ throw condition ([28]) or (j29|) is not satisfied 

if (not a: + ) < Oj or ("no* |^(i/,a;_) > 0J then 

|^ throw condition ([28f or ([29]) is not satisfied; 



26 end 
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Finally, in lines 15-25 we verify conditions (j2"51 - |2"9")) . Again we verify that for 
an element of grid v it holds G(min(z/), x+) > and G(min(p), xJ) < 0. This 
together with ^(v,x + ) > and $£(D,xJ) < proves (|2BH2SD- 

I 

As was mentioned earlier AlgorithmQ]is used to prove the existence of period 
two points curve for parameter values slightly above the first bifurcation - i.e. 
for parameters close to 02 G has three solutions close to one another. For these 
parameter values we found difficulties with verifying the existence of period-two 
points curve using C 1 -computations, only. 

For parameter values away from the bifurcation, where all eigenvalues of 
periodic orbit are well separated from the unit circle, we use Algorithm [2] based 
on the Newton interval method and the Krawczyk method, both requiring only 
C 1 -computations, and which verifies the existence of only one branch of period 
two points for P a . 

Before we present this algorithm we need to introduce some notations. Let 
fl = {(x, y, z) G M 3 : x = 0} be a Poincare section for ([3T)l and P a : fl-e-dl be 
corresponding Poincare map for a system with parameter value a. Notice that 
the trajectory can intersect n at a point (y, z) G n for which x' — —y — z is 
positive or negative (if it is equal to zero the Poincare map is not defined). Hence 
we have P% n = Pa, where P a is the Poincare map for section II = {(x,y, z) G 
R 3 I x — 0, x' > 0} and therefore period two points for P a correspond to period 
four points for P a . Let us define a map F a : LT 4 — > R 8 by 



(yi 






(2/2 


Z2) 


-P a (yi 


Z\) 


(2/2 


22) 




(2/3 




-Pa(V2 


Z2) 


(2/3 


Z3) 




(2/4 


z 4 ) 


~ Pa(V3 


z 3) 


.(2/4 


Zi)_ 




.(2/1 




~ Pa (2/4 


Zi)_ 



Algorithm [2] was used to verify the existence of a continuous branch of period 
two points for P a for a belonging to some interval. Sets Xi x Yi give the size of 
the neighborhood around a candidate periodic orbit on section II. Lines 4 to 8 
constitute a heuristic part and their task is to find a good candidate. 

Lemma 23 If Algorithm^ is called with its arguments [a*, a*], Yi x Zi, i — 
1,2,3,4 and g and does not throw an exception then there exists a continuous 
curve (yper, Zper) '■ [o*,o*] — > II such that (y per (a), z per {a)) is period two point 
for /',.. 

Proof: The existence of fixed point for P^ for all a G [a*, a*] is verified in lines 
12-16. Lines 17-18 guarantee that this is a period two point for P a , in fact a 
unique one in U. 

Uniqueness implies continuity on each a G Q and due to uniqueness and 
connectedness of the set B defined in line 19 we see that they agree on boundaries 
of a. 

I 

Proof of Theorem [21} 



31 



Algorithm 2: verification the existence of period two points branch. 
Data: [a*, a*] - an interval, Yi x Zj, i = 1, 2, 3, 4 - convex, compact sets, 
g - an integer 

Result: if algorithm stops and does not throw an exception then there 
exists a continuous branch of period two points for P a for 
parameter values a e [a», a*] 

begin 

Q < — .g-grid for [a*, a*]; 
foreach a G Q do 

a < — center (a); 

ui = (j/i, z\ ) < — find approximate period two point for P^ using 
standard Newton method; 

U2 = (v2,z 2 ) < — P a (yi,zi); 

U3 = (2/3,^3) < P_a{y2,Z 2 ); 

m 4 = (2/4,^4) < — P a (y3,z 3 ); 

C < — compute approximate value of DF a (ui, u 2 , u 3 , 114); 
if C is singular then 
|_ C < — Id ; 

£/ < — (ui + Y\ x Zx,u 2 + Y 2 x Z 2 ,u 3 + Y 3 x Z 3 ,u 4 + Y t x Z 4 ); 
u < — (ui,u 2 , w 3 ,u 4 ); 

#a = (fcl,a,fc2,a,^3,o,fc4,a) < 

IntervalKrawczykOperator(P a , C _1 , £7, w); 
if not if 5 C int U then 

|_ throw cannot verify the existence of period two point; 

if ki. a ("1 fc 3 , a 7^ then 

throw the unique fixed point for P| in fci !a is not necessary 
period two point for P a ; 

B < — Uaee« x 
if P is no£ connected then 

throw cannot verify if branch of fixed point curve is continuous 
on interval [a*, a*]; 

22 end 
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9i 


92 


53 


9x 


t 


(X, x Yi) • 10- 4 


[o 2 , 2.8325] 


6 


70 


600 


150 


0.88 


[-99.0216,97.95] x [-4,4] 


[2.8325,2.8326] 


10 


80 


300 


270 


0.7 


[-166.9,163.9] x [-4,4] 


[2.8326,2.8327] 


12 


60 


150 


320 


0.6 


[-214.59,209.62] x [-4,4] 


[2.8327,2.8328] 


17 


50 


120 


330 


0.5 


[-253.8,246.8] x [-8,8] 


[2.8328,2.8329] 


18 


50 


100 


350 


0.5 


[-287.757,279] x [-8,8] 



Table 4: Parameters of the Algorithm [TJ 



The existence of continuous curve of period two points on intervals (ai, 02] 
and [03,03] is a consequence of Theorem 1101 and Theorem I16[ respectively. Let 
a s = 2.8329. For parameter values [02, o a ] we verify the existence of period two 
points branch using Algorithm [T] and for parameter values a G [o a , GS3] we use 
Algorithm [2] 

We have ran Algorithm [1] five times with parameters listed in Table 2] (in 
each case the map is P a ). Since in each case the algorithm had stopped and did 
not throw an exception we conclude that in each interval of parameters listed in 
Table [3] there exist two continuous curves ci(a), and 02(a) of period two points. 
Sets Xi x Yi listed in Table [4] are chosen so that 

Y x Z c Xi x Yi C ••• C X b x Y 5 , (53) 

where Y x Z is the set used in the proof of the existence of first period doubling 
bifurcation. Observe also, that since we know that for holds P a2 (01(02)) = 
02(02) and 01(02), 02(02) G Y x Z C X\ x Y\ from Lemma [3] we obtain that 
{ci(a), 02(a)} is period two orbit for P a , for a <E [02, a s ], i.e. the whole interval 
of parameters covered by intervals listed in the first columns in Table [4j The 
uniqueness of period two orbit together with condition (|53|) implies that the 
curves are continuous on (oi, a s ). 

One can see that the total number of initial values for which we need compute 
third order derivatives of G,which is equal to the sum of g\g x over all rows in 
Table [4] is equal to 19350. The total time of computation of this step is ten 
hours on the Pentium IV, 3GHz processor. 

We have run Algorithm [2] with 74 different arguments listed in Table [5l 
We have chosen the parameters of the Algorithm [2] so that such 74 intervals 
[(a*)j, (a*)i\, i — 1) ■ ■ ■ , 74 cover the interval [a s , a 3 ]. Notice also, that for pa- 
rameters a closer to a s we need larger values of g since the hyperbolicity close to 
02 is very weak. The total number of subintervals used to cover interval [a s , 03] 
is 614450. In fact this is the longest part of the numerical proof. The total 
time of computation of this step is 53 hours on the Pentium IV, 3GHz proces- 
sor. Since in each case Algorithm [2] stops and does not throw an exception we 
conclude that on each subinterval [(o*)j, (o*)i] there exists continuous branch 
of period two points. We need to show that these curves glue continuously at 
a*'s. In fact, this algorithm returns an upper bound for this period two points 
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branch which is of the form 

74 

B=\jBi (54) 
»=i 

where B^s are defined in line 19 of Algorithm^ We verified that B is connected 
- this together with an information that for fixed a £ [0,3,0,3] there exists a 
unique period two point (y per (a), z per (a)) such that (a,y per (a), z per {a)) G B 
implies that the curve (y per (a), z per (a)) is continuous on [0^,03]. 

There remains to show the continuity of fixed point branch for parameter 
values a — a s and a = 03. For a = a s we know that there exist period two 
points ci(a a ), 02(0,3) which belongs to the last set listed in Tabled i.e. 

ci(o a ), c 2 (a s ) e W± = u + M ■ ([-287.757, 279] x [-8, 8]) • 1CT 4 

where uq and M define coordinate system close to first period doubling bifur- 
cation and are defined in (138ll3"9")l . On the other hand the estimation for period 
two point resulting from Krawczyk method used in Algorithm [2] is 

w 2 = (wlwi) 

W\ = [-4.7668051788293892557, -4.7667832743925968586] 
Wl = [0.052543190547910088861,0.052543238016254205369] 

One can verify that W2 C W\ which obviously means that a period two point 
(Uper(a>s), z per (a s )) £ W2 resulting from the Krawczyk method and Algorithm^ 
is one of the points ci(a s ), 02(0,3) resulting from Algorithm [T] Hence, the curve 
of period two points is continuous at a — a s . 

In a similar way we verified continuity at a — 03. From the Krawczyk method 
used in Algorithm [2] we know that (y per (aa), z per (a3)) is a unique period two 
point in the set 

w 3 = (w£,wi) 

W3 1 = [-4.5010116820607413146,-4.4996232549240025023] 
W3 2 = [0.043134233933640332703, 0.043140290681655812932] 

On the other hand from Theorem \W\ we know that for a — 03 period two point 
belongs to the set 

W± = u 2 + M 2 • (Y 2 x Z 2 ) 

where M2, u 2 , Y 2 , Z 2 define the set on which we verify the existence of second 
period doubling bifurcation. One can verify that W4 C W3 which proves that 
the branch of period two points is continuous at a — 03. 
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U = [-1, 1] • (1100, 3, 1000, 3000, 1100, 3, 1000, 3000) • 10~ s 


i 


[a*, a*] 


9 


Y x x Z x x Y 2 x Z 2 x Y 3 x Z 3 x Y 4 x Z 4 


1 


[2.8329,2.83291] 


14000 


U 


2 


[2.83291,2.83292] 


13500 


u 


3 


[2.83292,2.83293] 


12200 


3U 


4 


[2.83293,2.83294] 


11250 


3U 


5 


[2.83294,2.83295] 


10400 


3U 


6 


[2.83295,2.83296] 


9650 


W 


7 


[2.83296,2.83297] 


9050 


3U 


8 


[2.83297,2.83298] 


8500 


3U 


9 


[2.83298,2.83299] 


8000 


3U 


10 


[2.83299,2.833] 


7550 


3U 


11 


[2.833,2.83301] 


7200 


3U 


12 


[2.83301,2.83302] 


6800 


3U 


13 


[2.83302,2.83303] 


6500 


3U 


14 


[2.83303,2.83304] 


6200 


3U 


15 


[2.83304,2.83305] 


6000 


3U 


16 


[2.83305,2.83306] 


5700 


3U 


17 


[2.83306,2.83307] 


5500 


W 


18 


[2.83307,2.83308] 


5300 


3U 


19 


[2.83308,2.83309] 


5100 


3U 


20 


[2.83309,2.8331] 


4900 


3U 


21 


[2.8331,2.83311] 


4800 


3U 


22 


[2.83311,2.83312] 


4600 


3U 


23 


[2.83312,2.83313] 


4450 


3U 


24 


[2.83313,2.83314] 


4300 


3U 


25 


[2.83314,2.83315] 


4150 


3U 


26 


[2.83315,2.83316] 


4050 


3U 


27 


[2.83316,2.83317] 


3950 


3U 


28 


[2.83317,2.83318] 


3850 


3U 


29 


[2.83318,2.83319] 


3750 


W 


30 


[2.83319,2.8332] 


3650 


3U 


31 


[2.8332,2.8333] 


36000 


3U 


32 


[2.8333,2.8334] 


29000 


3U 


33 


[2.8334,2.8335] 


24000 


3U 


34 


[2.8335,2.8336] 


20000 


3U 


35 


[2.8336,2.8337] 


18000 


3U 


36 


[2.8337,2.8338] 


16000 


3U 


37 


[2.8338,2.8339] 


14000 


3U 


38 


[2.8339,2.834] 


13000 


3U 


Table\E\ continued on the next page 
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TaMe[5| continued from previous page 


i 


[a*, a*] 


9 


Yi x Zi x Y 2 x Z 2 x Y 3 x Z 3 x Y~ 4 x Z 4 


39 


[2.834,2.8345] 


59000 


3U 


40 


[2.8345,2.835] 


42000 


3U 


41 


[2.835,2.8355] 


17500 


15U 


42 


[2.8355,2.836] 


11000 


15U 


43 


[2.836,2.8365] 


8000 


15U 


44 


[2.8365,2.837] 


6200 


15U 


45 


[2.837,2.8372] 


2800 


30U 


46 


[2.8372,2.8375] 


3600 


30C7 


47 


[2.8375,2.838] 


5000 


30£7 


48 


[2.838,2.8385] 


3600 


30C7 


49 


[2.8385,2.839] 


2900 


30£7 


50 


[2.839,2.8395] 


2400 


30£7 


51 


[2.8395,2.84] 


2100 


30C7 


52 


[2.84,2.841] 


3700 


30C7 


53 


[2.841,2.842] 


3000 


30C7 


54 


[2.842,2.843] 


2500 


30£7 


55 


[2.843,2.844] 


2200 


30£7 


56 


[2.844,2.845] 


2000 


30C7 


57 


[2.845,2.846] 


1700 


30C7 


58 


[2.846,2.848] 


3100 


30J7 


59 


[2.848,2.85] 


2600 


30U 


60 


[2.85,2.86] 


11000 


30U 


61 


[2.86,2.87] 


6500 


30C7 


62 


[2.87,2.88] 


4700 


30J7 


63 


[2.88,2.89] 


3700 


30£7 


64 


[2.89,2.9] 


3100 


30C7 


65 


[2.9,2.95] 


4400 


150J7 


66 


[2.95,3] 


2500 


150C7 


67 


[3,3.1] 


3500 


150C/ 


68 


[3.1,3.2] 


2500 


150C/ 


69 


[3.2,3.3] 


2100 


150C7 


70 


[3.3,3.4] 


1800 


150C7 


71 


[3.4,3.5] 


1600 


150U 


72 


[3.5,3.6] 


1400 


150C7 


73 


[3.6,3.7] 


1500 


150J7 


74 


[3.7,03] 


2400 


150C7 



Table 5: Parameters of Algorithm [2l The initial set U is defined in the first line 
of the table. 
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7.1 Technical data. 



In order to compute Poincare maps P and P 2 with their partial derivatives 
we used the interval arithmetic [IE, MoJ, the set algebra and the C'-Lohncr 
algorithm (WZj developed at the Jagiellonian University by the CAPD group 
[CAPD]. The C++ source files of the program with an instruction how it should 
be compiled and run arc available at WIJ. 

All computations were performed with the Pentium IV, 3GHz processor and 
512MB RAM under Kubuntu Feisty Fawn linux with gcc-4.1.1 and MS Windows 
XP Professional with gcc-3.4.4. The computations took approximately three 
days. The main time-consuming part (over 63 hours) is the verification of the 
existence of connecting branch of period two points between first and second 
bifurcation. 
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